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Abstract 

In the past two decades, it has been estabhshed by high-resolution observations of early-type 
galaxies that their nuclear surface brightness and corresponding stellar mass densities are charac- 
terized by cusps. In this paper, we present a new spherical analytical model family describing mild 
cuspy centres. We study isotropic and anisotropic models of Osipkov-Merritt type. It is shown that 
the associated distribution functions and intrinsic velocity dispersions can be represented analyti- 
cally in a unified way in terms of hypergeometric series, allowing thus a straightforward comparison 
of these important global quantities for galaxies having underlying mass densities which may differ 
significantly in their degree of central cuspiness or radial falloff. 
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1 Introduction 

Since the early nineties of the last century, it has been established by observations of ground- and space- 
based telescopes that the nuclear surface brightness and corresponding stellar mass densities of early-type 
galaxies are characterized by cusps. The construction and study of galactic models incorporating cuspy 
centres has therefore been an active part of theoretical modeling, for which we can not give a full account 
here (see for instance [Tj, [S], [5], [1], [H], [H]). An important issue for its own has been the construction 
of analytical models describing cuspy densities, even more so since up to the eighties most of the ana- 
lytical models available included (flat) cores only. One of the first cuspy models, however, were given by 
[T5] and [12]. Further analytical spherical models with cusps have been presented later e.g. by [TU], [TO] . 
|21| and more recently by ^2 and 5 . However, many models in the literature which are able to capture 
varying degrees of cusps are rather inflexible with regard to the outer falloff behaviour of the density or 
vice versa. Additionally, the associated distribution functions must often be determined numerically, and 
analytical expressions, if they exist, are often limited to a few concrete values of the underlying model 
parameters. In this paper, we like to present a very general spherical model family which allows more 
flexibility in the central as well as in the outer radial shape of the associated mass density, starting from a 
family of non-singular, powerlaw-like potentials for the stellar component. All further intrinsic quantities 
can be calculated analytically, notably the distribution functions, for which an analytical representation 
for a large range not mere for a little subset of parameters is possible. This analyticity provides thus 
a straightforward study of the relation between the cuspiness of the density and the behaviour of the 
corresponding distribution function (and intrinsic velocity dispersion), than it would be without analyt- 
ical expressions at hand. To this aim, we consider an isotropic as well as an anisotropic parametrization 
of Osipkov-Merritt type for our models. The projected quantities follow straightforwardly, but must 
be determined numerically for our model family, unless the mass density has a (flat) core. The family 
presented below is able to model the cuspy centres of massive early-type galaxies and nucleated dwarf 
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elliptical galaxies. In addition, our family may also serve as useful input for numerical studies on the 
time-dependent evolution of galactic nuclei. 

This paper is organized as follows: In Section 2, we present our family of potentials and mass densi- 
ties. In Section 3, we deduce the distribution functions for isotropic and anisotropic Osipkov-Merritt 
parametrization. The intrinsic velocity dispersions are calculated in Section 4- In Section 5, we study the 
presence of a central supermassive black hole in our model and in Section 6 we present the conclusions. 
Appendix A contains some often used formulae. 



2 Model family 

We adopt the following family of spherical potentials 



- (1) 

with $0 and b positive numbers and <i>(r) ^ — as r ^ oo. In fact, this family comprises almost all 
non-singular powerlaw-like potentials found in the literature, most of which arc governed by one slope 
parameter. We choose the slope parameters according to < a < 2 and 7 > 0. In this paper, we like 
to study the self-consistent model originating from this potential. From Poisson's equation follows the 
corresponding mass density 

, , a76"'^$o (1 + + (1 - a7)r" 
p{r) - 



which is positive for < 1 and this restriction is imposed throughout. The cuspiness of the density is 
determined by the parameter a. Flat cores with a central density of p(0) = 37<I>o/27rG6^ are recovered for 
a = 2. Otherwise, there is a cusp with p{r) ^ 7^-2+0^ j^^ large radii, the density goes like p{r) ^ j.-'^-a-y ^ 
hence the degree of the outer falloff is governed by both parameters. The above family includes a lot 
of known models as special cases: For instance, the |T7] model is recovered by setting a = 2,7 = 1/2. 
For a — 1, the [12J model follows for 7=1. Other special cases obtained in the literature include 
a = 1/2, 7 = 2 and a = l,7 = /3 — 3 with /3 < 4, see [3T!- In order to recover the cusps of the models of 
[To] and [H], one had to put a — 1/n, n G N. On the other hand, the outer falloff is recovered by setting 
7 = 1/n. Both conditions at once can not be fulfilled to recover the full models of the above authors. 
However, we note that those models do not include mild cusps with < 2 — a < 1 whereas our density 
does. 

The associated cumulative mass function to ([T]) is given by 

. ^, , a7$nfc"'' 
M(r) = 



going for large radii as M(r) ~ hence only models with 07 — 1 have a finite total mass. In terms 

of the circular velocity 



..2 



this amounts to Vc ^ r'^"i/^ for large r. Thus, the circular velocity is Keplerian only for — 1, and 
decreases more slowly for a7 < 1. In the limit a7 — > it becomes constant, Vc — > const. However, for any 
fixed product G (0, 1), the increase in the cumulative mass is weaker as for the logarithmic potential 
([3]) or the isothermal sphere, where M(r) ^ r for large radii 0. Self-consistent models satisfying a7 < 1 
must therefore be cut off at some outer radius in order to provide a finite mass. On the other hand, due 
to its ability to reproduce constant or rising mass and velocity profiles at large radii, the potential in ([Ij 
may be also useful for modeling dark matter structures. In fact, as an example we refer to [30], where the 
intrinsic quantities for the stellar component were derived by assuming that a dark matter component 



^On the other hand, M(r) diverges only logarithmically for the Hubble-Reynolds or modified Hubble density profiles 
(see m). 



dominates the potential, the later of which is a special case of ([T]) with the parameters a 0, 7 = 1. 
In the forthcoming, it is advantageous to use dimensionless units: dividing ([T]) by — $0 and the mass 
density by a7(l + a)$o/47rG&^, we have for the (relative) potential 

= 77 ^ (3) 

and for the density (using the same notation) 

_ (1 + q)W + (1 - Q;7)r" 

^^^^ ~ 1 + a ^2-a(5a +^a)7+2 " ^ ^ 

These quantities will be used in subsequent calculations. We like to put our emphasis on intrinsic 
quantities which can be calculated analytically. Hence, the plots, which will be shown below, only display 
the distribution functions and intrinsic velocity dispersions, respectively. The corresponding projected 
quantities for the family in ([3]) & ([U must be determined numerically except in the case of cores with 
a = 2, where analytical expressions in terms of hypergeometric functions can be given as well. In any 
case, it can be shown easily that the surface brightness associated to (jl]) rises steeply with decreasing a 
as the projected radius tends to zero, since then there is more stellar mass concentrated in the nuclear 
region. This behaviour is more pronounced for larger values of 7. 



3 Distribution functions 



As is shown in this and the following section, the distribution functions (DFs) and intrinsic velocity 
dispersions for the above model family ([3]) & (|1]) can be calculated analytically and we are going to study 
their behaviour for varying cuspiness and outer falloff of the density Q . To this aim, we consider isotropic 
models with the DFs depending on the relative energy E — —E as well as anisotropic Osipkov-Merritt 
models, where they depend on £ and the angular momentum L via Q = £ — L^/ (2r^) (see [E] and [H]). 
The anisotropy radius is a free parameter and the anisotropy function for this parametrization behaves 
as (3{r) = r'^ /{r'^ + r^), hence the models are isotropic in the centres. 

The isotropic DFs f{£) for the model in ([S]) & (01) are calculated using Eddingtons's formula 
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by exploiting the fact that the density ^ can be expressed in terms of the potential ^ as 
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The function I{£) in ([5]) is then 
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The integrals in this expression can be determined analytically in terms of Beta functions B{a,b) = 
r(a)r(6)/r(a + 6) and hypergeometric series gFp(ai, 02, .., a^; 61, ^2, .., 6p; z) (see [TT] for definitions and 
properties) provided that 2/ a — n,l/"/ — m and 2/{nm) < 1 with n, m e N. For core models having 
a = 2, the only restriction on the value of 7, however, is to be < 1/2. We use now equation (fT7|) in the 
Appendix to calculate the integrals in which results into 
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where we defined for brevity 
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In order to calculate the derivative of ([7]), we use the general relation (fTS]) . Abbreviating 
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we finally arrive at the expression for the isotropic distribution function 
/(£) = (V87r2(l+a)7)-if-i/2 x 
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This functions involves powers of £ multiplied by hypergeometric series, the later of which may even 
reduce to simpler analytical functions of £ depending on the values for a and 7. The order of the 
hypergeometric functions is determined by 7. The DFs for models with cores, a = 2,7 < 1/2, simplify 
considerably and are given by 
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The function in fS]) is plotted in Fig. 1, first row, left plot, for the models a = 2,7 = 1/3; a = 
1,7 = 1/2; a = 2/3, j = I; a — 1/2,7 = 1/2- As a result of the finite depth of the central potential 
well, ^'(0) = £max — 1, the distribution functions diverge for £ ^ 1: As £ ^ 1, a steeper inner cusp 
corresponds to a stronger divergence in this limit because the system is then dominated by stars at small 
radii where the cusp dominates and this effect is therefore hardly affected by 7. On the other hand, the 
decrease of /(£) as £ ^ is larger for small values of 7. This is more pronounced if a is small as well 
because then the model is more centrally concentrated as a result of the cusp. 

Now we turn to the anisotropic models: The Osipkov-Merritt distribution functions for the model family 
in ([3]) & (HI) can be calculated accordingly from a similar relation as the one given in (O, (see f^), namely 
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using the auxiliary density Pa(r) 
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Figure 1: Log-log plots of the isotropic and anisotropic distribution functions and intrinsic velocity 
dispersions for the models a = 2, 7 = 1/3 (black solid), a = 1,7 = 1/2 (black dashed), a — 2/3,7 = 1 
(light grey), a — 1/2,7 — 1/2 (dark grey); first row: isotropic and anisotropic distribution function, 
second row: isotropic, radial and tangential velocity dispersion for r < 1, third row: isotropic, radial and 
tangential velocity dispersion for r > 1; G = 1 = 6, 7'a — 2 

The first term /(Q) is given by the expression for the isotropic DFs in ^ except that 8 has to be 
replaced everywhere by Q. The anisotropic DP ([9]) is plotted in Fig. 1, first row, right plot, using the 
same model parameters as before. For the same reason as in the isotropic case, the increase of fa{Q) for 
Q ^ 1 is dominated by the cusp parameter a. On the other hand, for fixed a the parameter 7 controls 
essentially the degree of the anisotropy in the sense that the model is more anisotropic for small values 
of 7. As a general result we see that the anisotropic DFs do not decrease as rapidly for Q ^ as do the 
isotropic DFs: It can be easily shown that the models approach the isotropic behaviour for large > 5, 
as expected. In contrast, for < 6 the anisotropic signature in fa{Q) dominates over a wider range in 
Q, whereas the increase for Q — )■ 1 remains quite unaffected. 



4 Intrinsic velocity dispersions 

The intrinsic velocity dispersions (VDs) for the isotropic models are derived from the usual relation 

r GM{s)p{s)ds 



as follows: Using equation (fT8|) for the radial range r > b after substituting s" = x, we get 
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whereas for the inner range r < b, suitable variable transformations and usage of equation ()16|1 results 
into 

(ha I ^Q\-7-l 
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In contrast to the DFs from above, the expressions here involve only special hypergeometric functions of 
the form 2Fi{a,b;c; z). In Fig. 1, second and third row, left plots, we show IV := y/ {vl{r)) for r <1 
and r > 1 respectively for the same model parameters as above. For r — > 0, the velocity dispersions even 
decrease for the weak cuspy models we are studying (this sounds counterintuitive, but that is a typical 
behaviour of such models as long as no additional central black hole potential is added, see e.g. [H]): 
For r ^ 0, the VDs converge to zero as r^" if a < 2. For a = 2, they are asymptotically constant, 
as expected: Note that the first hypergeometric function in (fTTj) then simplifies to (6^ + r^)/r^ and the 
second one to (6/r)^(l + (6/r)^)(l + 1/(2(7+ l)))? both denominators cancel with the factor in front. 
Moreover, the outer falloff also depends on the degree of the central cusp: The overall shape is flatter for 
increasing cuspiness, but this is already evident from the expression for the circular velocity given in 
It can be shown that the projected velocity dispersions exhibit the same overall behaviour with regard 
to the model parameters as do the intrinsic ones. 

Now we turn to the Osipkov-Merritt models: The intrinsic radial velocity dispersion is given by (see [7]) 

/ 2. ^^ 1 rl ( GM{s)p{s) , 1 w^. 

where the first integral was already evaluated in equations (jlOp and (fTTj) . The second integral, however, 
can be evaluated for r > 6 as 
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by using formulae ()18|) and (I16p after suitable substitutions, respectively. The intrinsic tangential velocity 
dispersion, on the other hand, is then simply given by 

for the respective radial range. In Fig. 1, second and third row, central and right plots, we show the 
radial and tangential velocity dispersion IV^ '■— \J {v'^{r))a and IVt :— {v\{r)) respectively, for the 
same parameters as for the isotropic models. Both dispersions decrease more slowly for small values of 7, 
i.e. for higher anisotropies. Concerning the overall shape, it can be shown that IV* falls off more rapidly 



for r > 6 than IV r- For increasing cuspiness, the shape of both velocity dispersions becomes flatter, 
although this behaviour is much less pronounced than it is for the isotropic models. The tangential 
velocity dispersion dominates over the radial velocity dispersion for > 5, whereas the opposite is true 
for Ta < b. Since {v'^{r)) 2{v'^{r))a for r ^ 0, the respective curves differ almost only by an overall 
factor of two for small radii. 



5 Models with central black hole 

In the presence of a central black hole of point mass Mbh, the density p(r) is not changed but the 
potential is modified according to 

**(r) = *(r) + ^^^i=:*(r) + ^. (12) 
9q r r 

The distribution fimction for the model ([2]), Jl]), (fT^ can be determined analytically only if the energy 
£ is large, i.e. close to the black hole. Since '^*{r) is no longer invertible with respect to r, one may 
perform the following variable transformation in I{£) of equation ([5]) (see e.g. |19j): 

US) = [' '^f'^^'^ = r^'^ "^^^^^ fi3i 

^ ' Jo rf*' V£~^ Jo du -^'(u)' ^ ^ 

where u = 1/r and u{£) is defined implicitly by ^!*(u(£)) = £. For large £ and small r (i.e. large u), we 
may approximate ^E"* pu and 

^ ^ (2 - a)b{ubr- + 2(l-^7)(l-«)^(^^)i-2.^ 
du 1 + a 

Inserting this into (fT5|) and using (fT7|) . the distribution function becomes then 
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which is valid for a < 1, but both terms are non- vanishing only for a < 3/4. However, both restrictions 
favour a cusp of 2 — a > 5/4 which is steeper as the r~^/^-cusp expected to be produced by the adiabatic 
growth of a black hole in the context of isotropic models. The distribution function is hardly affected by 
7 and the system is populated with more stars in the very centre, i.e. f'{£) is larger, with decreasing 
black hole mass fj,, as expected. 

In order to deduce the isotropic velocity dispersion in the presence of a black hole, we use the relation 

iv'^ir))' = {vUr)) + ^J^^ds, 
p{r) 

where the first term is the velocity dispersion from cqu. (jlOp and (jlip . and the second term can be 
evaluated as 
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for r > 6 by using (fT^ and respectively. In Fig. 2, we show IVj, := \J {v'^{r))* for r < 1 and r > 1. 
Now the velocity dispersions rise steeply for r ^ in contrast to the previous case without black hole. It 
can be also shown that this rise at small radii is more pronounced if /i is increased but a is fixed together 
with a more slowly falloff in the outer parts. The same overall behaviour is found for the corresponding 
projected velocity dispersions. 



Figure 2: Log-log plots of the intrinsic velocity dispersion for models with black hole: a = 2/3,7 — 1 
(light grey), a = 1/2,7 = 1/2 (dark grey), ^J^ = 0.001, G = 1 = 6 



6 Conclusions 

In this paper, we considered a family of non-singular potentials falling off as 1/r or more slowly at large 
radii. The associated self-consistent mass density incorporates flat or cuspy nuclear regions together with 
a flexible falloff behaviour at large distances from the centre. The corresponding distribution functions 
and intrinsic velocity dispersions can be represented analytically in terms of hypergeometric functions. 
This allows a straightforward comparison between models for galaxies having different central and outer 
shapes in the mass density. We restricted ourselves to isotropic and anisotropic models of Osipkov-Merritt 
type. It is shown that the anisotropy affects the distribution functions only outside the central parts where 
they do not fall off as rapidly as the isotropic ones, whereas the increase for large arguments is dominated 
by the cusp parameter in both cases. Moreover, the velocity dispersions decrease more rapidly for the 
less anisotropic models and their shape is flatter for increasing cuspiness. The presence of a central point 
mass potential, mimicking a massive black hole, is also studied. It is shown that the velocity dispersions 
rise steeply at small radii for increasing black hole mass, which in the same time leads to higher values 
of the velocity dispersion over a wider radial range. 
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A Formulae 

The formula for the derivative of the general hypergeometric series is applied in Section 5, 

pFq{ai,a2, ...,ap;bi,b2, ■■■,bq;z) = "^^^"""^ pFq{ai + l,a2 + l, flp + l; 6i + l, 62 + I, ...,6, + l;z). (15) 

UZ 0lU2---0q 

In Section 4, the transformation formula for the special hypergeometric function is used 

2Fi{a,b;c;z) - (1 - z)"" 2F1 (^a,c - b;c; . (16) 
Following integral relations are used in the text, see |Tl] : 

x^^^u - x)^'-^{x"' + (i"' fdx = |3"'^u^'+''-^B{^l, v) x 

/ t^J^+l v + m—l fi + v + fi + v + m — 1 / u \ ™\ 

m+i^m — A, — , ; , '^^[15] 

\ m m m m m m \p J J 

if Re{n) > 0, Re(t/) > 0. 

x-^{x + (3)''{x - u)''-^dx u^+^-^BiX 2F1 X - - i^; X - 1^; (18) 

if < Re(/i) < Re{X-iy). 







